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ABSTRACT 

We use data from the literature to constrain theoretical models of galaxy formation. 
We show how to calculate the dimensionless spin parameter A of the halos of disk 
galaxies and we compare the distribution of A with that observed in cosmological 
iV-body simulations. The agreement is excellent, which provides strong support for 
the hierarchical picture of galaxy formation. Assuming only that the radial surface 
density distribution of disks is exponential, we estimate crudely the maximum-disk 
mass-to-light ratio in the /-band and obtain (Tj) SS 3.56/i, for a Hubble constant of 
100/ikm 8 _1 Mpc _1 . We discuss this result and its limitations in relation to other inde- 
pendent determinations of Tj. We also define a dimensionless form of the Tully-Fisher 
relation, and use it to derive a value of the baryon fraction in disk galaxies. For galaxies 
with circular velocity v m > 100km s _1 , the median value is = 0.086( T//3.56/i). 
Assuming that the gas fraction in galactic halos is at most as large as that in clusters, 
we also conclude that (T/) < 2.48/i _1/2 . 
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1 INTRODUCTION 

The fact that the rotation curves of spiral galaxies are rather 
fiat is usually taken to imply the presence of an extended 
halo of dark matter (e.g. Freeman 1970, Persic & Salucci 
1991). If such dark matter exists, then hierarchical mod- 
els of galaxy formation (White & Rees 1978) are a natural 
consequence of gravitational instability. In these models the 
standard picture of disk formation is that gas, which is ini- 
tially distributed in the same way as the dark matter, cools 
and settles into rotationally supported disks at the centres 
of dark matter halos (Fall & Estafthiou 1980). 

Dalcanton, Spergel & Summers (1997), and Mo, Mao & 
White (1998) (MMW) have shown recently that such a pic- 
ture can reproduce some of the broad properties of observed 
disk galaxies. In particular they show that the Tully-Fisher 
relation (Tully & Fisher 1977) can be understood very sim- 
ply. However, as was pointed out by Courteau & Rix (1997), 
the Tully-Fisher relation is incompatible with the maximum- 
disk hypothesis (Carignan & Freeman 1985) if disk galaxies 
have universal mass-to-light ratios (de Jong 1996). 

Great interest in the Tully-Fisher relation as a distance 
indicator (Giovanelli etal. 1997 and references therein) has 
led to many observations of disk galaxies, and large samples 
are now available (e.g. Mathewson & Ford 1996; Courteau 
1996, 1997). In this paper, we examine the observational 
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constraints which can be placed on the standard picture of 
disk formation. 

In the next section we review the properties of exponen- 
tial disks and the scenario of disk formation. We start from 
the standard assumption that the disk mass-to-light ratio T 
is universal (or at least does not vary strongly with surface 
brightness), and calculate various components from observ- 
able quantities. In Sectionta we describe the observations of 
disk galaxies. In Section klwe show the results of applying 
the theory of Section [| to the data. In Section ^| we discuss 
the implication s of our results and draw conclusions. In par- 
ticular, Section 5.1 is devoted to a discussion of independent 



determinations of T. 



2 DISK FORMATION 



2.1 Exponential disks 



The luminous disks of spiral galaxies are commonly mod- 
elled by an exponential surface brightness distribution: 



H(R) 



2-kRI 



exp(-i?/7? d ) 



(1) 



where R is the usual cylindrical radius, R c \ is the exponen- 
tial scalelength, and Ld is the total luminosity of the disk. 
Here we collect some notation and a number of useful results 
relating to exponential disks. 

The disk has a mass M<j , and a mass-to-light ratio T in 
solar units. Thus the surface mass density of the disk is 
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E(i?) = n(R)T = 



exp(-J?/i? d ). 



(2) 



The gravitational potential in the disk $ is conveniently 
decomposed into contributions from the disk and a halo: 



(3) 



(We use the subscripts 'd' for 'disk', and 'h' for 'halo' 
throughout.) We assume for the present purposes that the 
halo is spherical, and usually we think of it as being com- 
posed of dark matter, but it may also contain a stellar com- 
ponent (e.g. the 'bulge' of an earlier type spiral). 

The speed of test particles on circular orbits v c as a 
function of R is given by 



v!(R) 



-R 



9£ 
dR' 



(4) 



We shall refer to v c (R) as the rotation curve of the system. 
The rotation curve as measured in HI (apart from small 
contributions from turbulent motion) is thought to be a good 
measure of the true rotation curve as long as the system is 
axisymmetric. 

We characterise the self gravity of the disc through the 
the dimensionless quantity 

(5) 



where Vm is the maximum value of v c , and G is the gravita- 
tional constant. The rotation curve of an isolated exponen- 
tial disk ($h = 0) is given by Freeman (1970). An isolated 
disk has e m = £d ~ 0.63 and a disk embedded in a halo has 

The directly observable quantity corresponding to e m is 



ez = 



(GL d /R d y/i ' 



which is related to e m by 
ei = e m T 



(6) 



(7) 



The quantity ej 3 has the units of a mass-to-light ratio, and 
indeed it is a measure of the total mass (including dark mat- 
ter halo) contributing to the rotation curve. Let us define a 
quantity 



T tot (7?) 



v 2 c(R)R 
GL(R) ' 



(8) 



which measures the total mass-to-light ratio as a function 
of radius. For an isolated disk it is a constant (= T) and 
with an extended dark halo it increases with radius. The 
maximum rotation velocity in general occurs at R = _R max > 
Rd, and the luminosity enclosed is L max < id, hence ef < 
T tot (7? max ). For an isolated disk i? max = 2.2Rd and I/ max = 
0.65L d , so T = T tot (i? max ) = 3.4e?. 



2.2 Disk Formation Model 

Here we reproduce a simple model where the dark matter 
halo is assumed to be a singular isothermal sphere and disc 
self-gravity is neglected (cf MMW). In this model, the disk 
scale length is 



Rd — — =Ai?h 

\/2 



(9) 



where A is the dimensionless spin parameter, R^ is the virial 
radius of the halo. The halo properties are given by 



M h = 



«h-Rh 



Rh = XjT 



(10) 



where Mh is the mass and Vh the circular velocity of the halo, 
Ho — 100/ikm s _1 Mpc _1 is the Hubble constant and x i s a 
dimensionless constant. For discs assembled at redshift z, 
MMW argue that x ~ 0.1H /H(z), where H(z) is the Hub- 
ble constant at z. We can treat x as an adjustable parameter 
and derive its value from observational data. Whenever we 
need a n umerical value for \ we use the one derived in Sec- 
tion 4.2 (x = 0.049). The maximum rotation velocity of the 
disk is v m » v h . 

When more realistic halo profiles are used and disc self- 
gravity is taken into account, the relation between Rd and 
Vh is slightly modified from that given by equations (|^) and 
( [i"c)| ) (see MMW for details) . The constant part of such mod- 
ifications can be taken into account by the constant x- For 
simplicity, we will ignore all high order effects. 

The disk central surface density is given by 

where is the fraction of halo mass that settles into the 
disk. We now define a dimensionless ratio of observables 

Ho v m 1 



\ z — 



h irGfjio T 



1.47 x 10" 



( *2 

V 200km s~V 



lOOL pc" 
Mo 



(12) 



where fio = Eo/T is the central surface brightness. Defining 
also the structural parameter 



rridh 1 
~ X 



(13) 



we combine equations (H-|13|) to obtain an expression for A: 



A 



From equations (^) and ([nj) we find that 

_L A_ 

V2m d 



(14) 



(15) 



which is a dimensionless description of disk formation: the 
smaller the spin parameter A, the more concentrated the 
disk, and the more self-gravitating it is (smaller e m ). We 
can write equation ( |l5| ) in terms of observables as 

v2x Om' 

Equation (|l(j) is just the dimensionless form of the Tully- 
Fisher relation. To see this we write the Tully-Fisher relation 

as 

H 2 L = Avl, (17) 
where the value of A can be derived from equations (j{^-|ll]): 

(18) 



A 2 Hq 

A = a ™ X hG 



Equations jl^ ) and (^) both relate observable quantities via 
the same constant of proportionality (i.e. a m x 2 ), and hence 
( [fif ) is a dimensionless form of the Tully-Fisher relation ( |l7|) . 
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Figure 1. The relationship between A m and the central surface 
brightness fio (in the I band) of the sample of spiral galaxies of 
Mathewson & Ford (1996). Barred galaxies are shown as triangles. 



10.0 - 




Figure 2. The dimensionless Tully-Fisher relation (equation [l6[) 
for the MF data set. On this plot, a fixed value of e m corresponds 
to a horizontal line with amplitude £^(Yj/i -1 ). The horizontal 
solid line shows the value of e m for a self-gravitating disk (e m = 
0.63) with Yr = 3.56h. The two dashed lines contain 90% of the 
data with slope and scatter predicted by the Tully-Fisher relation, 
within the framework of MMW (see text). Barred galaxies are 
shown as triangles. 



An important parameter in the above model is the di- 
mensionless spin A of the system before disk formation. In 
the hierarchical model of galaxy formation, much is known 
about its expected value and distribution (Cole & Lacey 
1996, Warren etal. 1992, Lemson & Kauffmann 1998). We 
can use equations (^) and (^) to write 



V2 X 



(19) 



the right hand side of which contains only observable quan- 
tities (and the constant x). Thus, this equation can be used 
to derive the distribution of A from observations once the 



value of x is fixed (see Section 1.2) 



3 OBSERVATIONS 

To compare theory with observations we use the dataset of 
Mathewson & Ford (1996) (MF) which has rotation veloc- 
ities and /-band photometry for a sample of nearly 2500 
Southern spiral galaxies selected randomly from the ESO- 
Uppsala catalogue (Lauberts 1982). The majority are rela- 
tively late types: of those 2275 for which Hubble types are 
given, 1055 are Sc or Sbc; 814 are Sb; and 5 are Sa. We 
convert the published photometric quantities to Rd and the 
central surface brightness fio by assuming an exponential 
profile. Details are given in an Appendix. 

The observed Tully-Fisher relation (Giovanelli et al. 
1997; Shanks 1997): 

Mi-5logh = -(21.0Q±0.02)-(7.68±Q.13)(log W-2.5), (20) 

where W is the inclination-corrected width of the HI line 
profile. The MF data has a Tull-Fisher relation which is 
consistent with this, albeit with larger scatter than in more 
carefully selected samples. 

The maximum of the rotation curve v m is given by 
Vm = W/2. Figure ^ shows the relationship between A m and 
the central surface brightness fio of the MF data set. The 
correlation is expected, since halos with smaller A should 
form more compact disks with higher Eo. This figure illus- 
trates that A and fj,o are almost interchangeable, despite the 
fact that in principle the correlation could have been washed 
out by scatter in v m . The tight correlation between A m and 
fio therefore implies that the disk central surface bright- 
ness is determined mainly by halo spin parameter rather 
than by halo circular velocity, as is expected in the disk 
formation model (see §2.2). From equations (|^)-(fi"l|) we see 
that the disk central surface density scales with Vh and A as 
Eo oc «h/A 2 . The dynamical range of for disk galaxies is a 
factor of about 3, while that for A 2 predicted for dark halos 
is a factor of about 15. 

In Figure ^, we show e; as a function of A m for the MF 
data. The extra factor of h in the abscissa makes the plotted 
quantities independent of the Hubble constant. The figure 
reveals a marked correlation between e; and A m which corre- 
sponds to the dimensionless Tully-Fisher relation (equation 
|li|). The slope is unity provided a m is independent of A m . 
The figure shows that the slope is consistent with unity for 
the majority of galaxies. At the large A m end the observed 
ei is slightly higher than the model prediction. There ap- 
pears to be a break in the slope of €i versus A m at A m ~ 0.1. 
According to the model in Section 2.1 a higher value of e; 



would mean that the combination rridh/T is lower for low- 
surface-brightness galaxies. This would result from a lower 
star formation efficiency (hence higher T), as suggested by 
observations (e.g. McGaugh & de Blok 1997). 

For a given A m , the scatter in e; is determined by the 
scatter in the Tully-Fisher amplitude A. In Figu re B l we 
overlay the predicted slope (i.e. unity, see equation |l6[) and 
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Figure 3. The distribution of T/ for the MF data, on the as- 
sumption that all disks are exponential and self-gravitating. The 
solid line is the whole sample, and the dashed line is for the barred 
subsample. The error bars are estimated using a standard boot- 
strap method. 



scatter on top of the data points with the normalisation 
chosen to reproduce the observed median value of e; in the 
range A m £ (.05, .06). The lines are derived from the 5% and 
95% quartiles of A of the data — the equivalent scatter in A 
is a factor of ss 2.4 or 0.95 magnitudes. The predicted slope 
and scatter are consistent with the observed data points. 
Thus, the data are consistent with the assumption that the 
value of a m is independent of A m for the majority of these 
galaxies. This is an important result, because it means that 
the total mass-to-light ratio, Mh/Ld = T/mj oc a m , is a 
constant for the majority of galaxies. We will return to this 
point in the next section. 



4 CONSTRAINING GALAXY FORMATION 

In this Section, we use the observational data to constrain 
the disk model described in Section [T^. We concentrate on 
the three most important parameters in the disk model: the 
mass-to-light ratios of disks, the spin parameters of halos, 
and the baryon fractions in disks. 



4.1 Mass-to-light ratios 

Assuming that the disks are exponential sets a limit on Tj: 
they should all have e m > 0.63, the value for an isolated 
disk. Thus for each galaxy Tj Si (ej/0.63) 2 , which is a crude 
estimate of the maximum-disk mass-to-light ratio. Figure 
|§] shows the distribution of Tj calculated in this way for 
the MF data, and for the barred subset. In Figure @ the 
solid horizontal line marks the median of e; in the data of 
Mathewson & Ford (1996). This corresponds to e m = 0.63 
for T/ = 3.56h. 

A conservative upper limit for the average disk galaxy 
would be (Tj) < 3.56/i. This limit is also consistent with 
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Figure 4. The distribution of A for the MF data (solid line with 
error bars). The best fit to a log-normal distribution (A = 0.05, 
a = 0.36) is shown as a dotted line. The error bars are estimated 
using a standard bootstrap method. 



independent measurements of Tj (see Section 5.1 for de- 
tails). However, since e; is higher for low surface brightness 
galaxies, a higher T/ is still allowed for these galaxies with- 
out violating the constraint e m > 0.63. Similarly, for high 
surface brightness galaxies, a lower Tj is required to avoid 
violating the constraint e m > 0.63. Note that if Tj were not 
universal, then the Tully-Fisher relation would require that 
md a Tj. In the opposite case, if Ti were universal, a more 
stringent upper limit on Tj would be required, in order to 
accommodate high surface brightness galaxies. 

4.2 Spin parameters 

The distribution of A can be determined from JV-body sim- 
ulations of hierarchical clustering, and is found to be log- 
normal with mean A « 0.05 and standard-deviation a ~ 0.5 
almost independent of cosmology (Lemson & Kauffmann 
1998, Warren etal. 1992). If observations of disk galaxies 
are a fair sample of dark matter halos, then the distribu- 
tion of A m should be closely related to the distribution of 
A. The distribution of A in the MF data (calculated from 
equation ^) is shown in Figure ^. The value of A has been 
fixe d at 0.05 by choosing \ = 0.049. As discussed in Section 



2.2 



X ~ 0.1-ffo/ H(z) for isothermal halos. Using more real- 
istic halo profiles and taking into account disk self gravity 
reduces the value of i?d for a given t> m , which corresponds 
to a reduction of the value of X- This value of x found here 
is actually in good agreement with the detailed modelling of 
MMW. 

The distribution is remarkably close to log-normal with 
a = 0.36 ±0.01. The value of a derived from the data is 
smaller than the value, a = 0.5, given by N-body simula- 
tions of dark halos. The reason for this discrepancy may 
be due to the fact that the observational sample is biased 
against both low-surface brightness galaxies, which are as- 
sociated with high-spin systems according to the disk model 
considered here, and early type spirals, which are associ- 
ated with low-spin systems. The main effect of the selection 
function is to rule out galaxies below some threshold in sur- 
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4.3 Baryon fractions 



Figure 5. Shows mj = \a m T /h versus e; for the MF data, as- 
suming that Xi is given by the crude upper limit Tj = (e;/0.63) 2 . 
Barred galaxies are again plotted as triangles. 




Figure 6. The MF galaxies which have B—V colours in the 
ESO/Uppsala catalogue. There is a moderate scatter in B—V 
but no trend with A m . 



face brightness (corresponding to 23.5mag/arcsec 2 ). There 
is also a distance dependent component, selecting galaxies 
with a range of sizes which depends on their distance. To 
check that the selection function does not severely affect 
our results we constructed Monte-Carlo samples using the 
model in Section 2.1 and applied the selection function to 
them. The selection process reduced the apparent width of 
the A distribution; the best fit to the data came from a dis- 
tribution with a\ ~ 0.4. The values of \ an d n-d (see below) 
are derived from A and are little affected. One feature of the 
surface brightness threshold is that it tends to cut off the 
low-luminosity (low-t> m ) galaxies in the Tully-Fisher rela- 
tion, reducing the slope from the theoretical value (3) closer 
to that observed (« 2.5 in the MF sample). 



Once x is fixed 
ing equation (P 
100km s _1 we find that 



a m can be estimated for each galaxy us- 
. Restricting to those galaxies with v m > 



<a m ) = 0.47 . (21) 

From equations © and @ we have the total mass- 
to-light ratio in the /-band: 

1 



{M h /L I ) 



= 43ft. 



(22) 



The observed mean luminosity of the Universe derived from 
recent redshift surveys of galaxies gives a critical mass-to- 
light ratio {M/L) w 1500/i(M/L)© to close the universe 
(see e.g. Lin et al. 1996). Thus the mass associated with 
individual galactic halos gives fi ga i = 43/i/1500fr ~ 0.03, 
which is a small fraction of the total mass in the universe. 
Using equations (jl3|) and ( pl| ) 

<™ d > = - 086 ^- w 

The result is shown as a scatter plot in Figure [| using the 
crude upper limit T/ = (q/0.63) 2 . The gas fraction ob- 
tained for X-ray clusters is / gas ~ 0.06ft _3//2 (e.g. Evrard 
1997). This fraction is usually considered to be equal to the 
mean value for the whole Universe (e.g. White etal. 1993), 
and therefore should be at least as large as that in galaxies. 
Equation ( p3| ) together with the gas fraction in clusters thus 
imply that 



(Ti) <, 2.48 h~ 1/2 . 



(24) 



which is an new upper limit independent of that derived in 
Section [0]. Clearly, the observational constraint on the frac- 
tion of baryons in galactic halos has important implications 
for the formation of galactic disks. 



5 DISCUSSION 



5.1 Independent determinations of T 



The disk mass-to- light ratios we derive in Section 4.1 are 
formally upper limits, and based only on rotation curve data. 
Thus they are consistent with mass-to-light ratios derived by 
other authors using the maximum disk hypothesis (Carignan 
& Freeman 1985, Palunas & Williams 1998). 

Limits on the mass-to-light ratio of the Galactic disk 
in the solar neighbourhood can be derived from a combi- 
nation of kinematic measurements and star counts. Kuijken 
& Gilmore (1989) derive a local surface mass density in the 
disk of 40Mq /pc 2 , and star counts give a V^-band luminosity 
density of 15L /pc 2 (Gould, Bahcall & Flynn 1996). Divid- 
ing mass by light we obtain Tv = 2.67, and thus Ti ~ 1.9 
(assuming V — I = 1.0). This number is independent of h, 
but is comparable with the median max imum-disk value for 
the MF galaxies derived in Section 4.1 for h <; 0.5. 



Mass-to-light ratios can also be derived from pure stel- 
lar population synthesis arguments, although there is al- 
ways some uncertainty arising from the poorly known initial 
mass function (IMF) , particularly from the low-mass cut-off 
in the IMF. A few stellar population models are available 
(e.g., Bertelli et al 1994; Bruzual & Chariot 1993; Worthey 
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1994). For a Salpeter IMF, the mass-to-light ratio derived 
by various authors appear to agree within an accuracy of 
25% (Chariot, Worthey & Bressan 1996). The predicted Tj 
depends on the metallicity and age of the stellar population. 
For a stellar population with age between 5-12 Gyr and with 
a solar metallicity, Ti is between 0.9-1.8 for a constant star 
formation rate (cf Table 3 in de Jong 1996). For an exponen- 
tial star formation law, the mass-to-light ratio is about 20% 
higher. The predicted values are in good agreement with the 
values derived in this paper. In the comparison, we have ne- 
glected the uncertainty due to dust, since the Tully-Fisher 
studies already attempt to correct for its effect. It has been 
argued (de Jong 1996) that dust reddening probably plays 
a minor role in the colour gradients in disk galaxies. Nev- 
ertheless, the dust correction remains a nuisance in these 
comparisons. 

The mass-to-light ratio of extragalactic disks can also 
be measured directly from kinematic studies such as that 
described by Bottema (1993). The measurement of T relies 
on the relation between vertical velocity dispersion <r z , sur- 
face density E and vertical scale height zo- Essentially, the 
larger the value of E, the hotter a disk has to be at constant 
zq. Using a small sample of very bright galaxies, Bottema 
(1997) derives a value of Ti = (1.7 ± 0.5)h, which is con- 
sistent with the upper limits we derived in Sections [O] and 



4.3. It is important to extend the range of data analysed 
by Bottema (1997). This requires a good HI rotation curve 
for each galaxy, and high quality spectroscopy at least along 
the major axis of the galaxy. Ideally one would choose the 
brightest members of a large pre-defined sample (such as 
that of Mathewson & Ford 1996) and follow them up with 
a high spatial resolution spectrograph. Kinematic informa- 
tion in more than one dimension is also advantageous since 
it removes some of the uncertainties in the deprojection of 
the velocity ellipsoid. A number of two dimensional spec- 
trographs are due to come on line shortly, and these may 
be well suited to the problem. The biggest uncertainty in 
the determination of Ti will remain that associated with 
the value of zq. Efforts should therefore be made to anal- 
yse as many edge on galaxies as possible to try to improve 
on existing determinations of the distribution of zo / Ra (e.g. 
Barteldrees & Dettmar 1994). 



5.2 Uncertainties in the results 

The constant x m equation (^) might in principle have been 
different owing to angular momentum loss. The fact that the 
derived value is close to the one expected in the disk forma- 
tion model implies that the gas should not have lost much 
angular momentum during disk formation. This conclusion 
was also reached by Fall & Estafthiou (1980) and MMW 
based on the observed disk scale lengths of local galaxies. 
We confirmed this using our Monte-Carlo simulations: de- 
creasing the value of \ led to disks with smaller 7?d and 
larger /i at fixed w m compared with the MF sample (L<j is 
fixed by the Tully-Fisher relation). 

According to MMW, when realistic halo profiles are 
used and disk self- gravity included, the factor of l/y/2 in 
equation (^) and the constant \ m equation (|lfj|) should in 
principle be replaced by quantities that depend on the spin 
parameter and concentration of galactic halos. The details 
of these dependences are uncertain, because they require ac- 



curate knowledge on halo density profiles, gas settling pro- 
cesses and star formation feedback. To derive observational 
constraints on these details of disk formation requires de- 
tailed decompositions of individual galaxies into different 
components as well as accurate measurements of disk rota- 
tion curves. Such analysis is not possible with the data set 
used here. 

the limit on the disk mass- 



As discussed in Section 4.1 



to-light ratio was derived assuming that Ti is independent 
of no (or of A m , see Figure [j]). This assumption is, as we ar- 
gued, consistent with the observational data and with inde- 
pendent measurements of Tj (cf. Figure []). Since the Tully- 
Fisher amplitude A derived from the data is quite indepen- 
dent of fio for the majority of galaxies (cf. Figure ^f, any 
trend in T/ with /io has to be compensated by a similar 
trend in m&. This would happen if disks with lower surface 
brightness contain larger amount of gas but have a lower star 
formation efficiency. Such a trend is suggested by the obser- 
vations of McGaugh & de Blok (1997) where gas fraction in 
low surface-brightness disks is compared with that in high 
surface-brightness disks. However, for most of the galaxies in 
the MF sample, disk masses are expected to be dominated 
by stars rather by gas, and so the trend in the gas fraction 
with surface brightness should not induce a significant trend 
in the disk mass-to-light ratio. 



5.3 Barred galaxies 

It is well known that isolated disks are violently unstable. 
Since disks in halos with lower spin parameters are more 
compact and more self-gravitating, they are more prone to 
global instabilities. One obvious possibility is that globally 
unstable disks turn into barred galaxies. If so one might 
expect that barred galaxies to have systematically smaller 
A m and larger /io- No such trend is seen in Figure ^ where 
barred galaxies (triangles) seem to be randomly drawn from 
the galaxy population. Two possibilities occur to us. First, if 
T is universal, then global instabilities in disk galaxies must 
be switched off by some mechanism which is independent of 
Ho. Such a mechanism might involve central mass concen- 
trations (Toomre 1981) which stabilise a disk by interfering 
with transmission of density waves through the centre (Sell- 
wood & Moore 1998). Second, if T depends on fio, then it 
must do so in such a way as to make all galaxies equally 
susceptible to bar formation. Bars could then be formed by 
interactions between galaxies (Noguchi 1996), or between 
galaxies and their dark matter halos. (N.B. the Tully-Fisher 
relation then requires that m,d depends on no-) 
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M = exp(q) - (a + 1) , 2g v 
/i a 2 

Given fi/fj, we can solve equation (|2^) numerically to find a. 
Then we have Ra = R/a, and /io = /iexp(a). 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science I^TpX style file. 



APPENDIX 

Here we describe how values of Rd, /io and L were de- 
rived from the published quantities of Mathewson & Ford 
(1996). The published data list total magnitudes, and face- 
on corrected isophotal quantities: average surface bright- 
ness (jl) and isophotal diameter (R) at /i corresponding to 
23.5mag/arcsec 2 . Assuming an exponential disk we have a 
surface brightness profile 

R 

Li(R) = fioexp(-a), a = — , (25) 
and 

L(R) = L d (l- exp(-a) (a + 1)) . (26) 

From equation (^) the average surface brightness inside ra- 
dius R is 

_ = i^ = 2^ 
itR* 1 a z 

where we have used Ld = 2-k[iqR\. Combining equations 
( pif ) and (p7|) we obtain 
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